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Abstract 

We report the results of a study of the exclusive charmless semileptonic ■k~£^i' decay 

undertaken with approximately 227 million BB pairs collected at the T{4S) resonance with the 
BABAR detector. The analysis uses events in which the signal B mesons are reconstructed with a 
novel loose neutrino reconstruction technique. We obtain partial branching fractions in 12 bins of 
q^, the invariant mass squared, from which we extract the f^{q'^) form factor shape and the 
total branching fraction: B{B'^ ■k~1'^v) = (1.44 it O.OSstat ± O-lO^ysj) x 10~^. Based on a recent 
theoretical calculation of the form factor, we find the magnitude of the CKM matrix element \Vub\ 
to be ^4.1 lb Q.2stat i ^■^gysU-^A'F'p) ^ 10^^, where the last uncertainty is due to the normalization 
of the form factor. 
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1 Introduction 



The precise measurement of \Vub\, the smallest element of the CKM matrix will strongly con- 
strain the description of weak interactions and CP violation in the Standard Model. 

The measurement of \ Vub\ requires the study of a 6 — > n transition in a well-understood context. 
Semileptonic b ulv decays (here, I stands for an electron or a muon) are best for that purpose 
since they are much easier to understand theoretically than hadronic decays, and they are much 
easier to study experimentally than the less abundant purely leptonic decays. 

For — > ■K~i'^v decays,^ the theoretical description of the quarks' strong interactions is 
parametrized by a single form factor, f~^{q^), where q'^ is the squared invariant mass of the i^v 
system. Only the shape of f^{q'^) can be measured experimentally. Its normalization is provided 
by theoretical calculations which currently suffer from relatively large uncertainties and, often, do 
not agree with each other. As a result, the normalization of the f'^{q^) form factor is the largest 
source of uncertainty in the extraction of \Vuh\ from the 7r~£^u branching fraction. Values 

of f^{q^) for B^ ir^i^v decays are provided by unquenched 0131 and quenched |3| lattice QCD 
(LQCD) calculations, presently reliable only at high q^ (> 16 GeV^/c^), and by Light Cone Sum 
Rules calculations j6j (LCSR), based on approximations only valid at low q'^ (< 16 GeV^/c^), as 
well as by a quark model [Z]. The QCD theoretical predictions are at present more precise for 
B^ TT^l^u decays than for other exclusive B Xulv decays. Experimental data can be used 
to discriminate between the various calculations by measuring the f~^{(f') shape precisely, thereby 
leading to a smaller theoretical uncertainty on | V^b \ ■ 

The present analysis of B^ 7r~i~^i' decays aims to obtain an accurate measurement of the 
f~^{q'^) shape in order to extract a more precise value of \Vub\ from the measurement of the total 
B^ — > 7r~i~^i' branching fraction, B{B^ — > 'k~1'^i'). To do so, we extract the B^ — > t:~1'^v yields in 
12 bins of q^ using a loose neutrino reconstruction technique and g^-dependent cuts. The quantity 
q^ denotes the uncorrected measured value of q^ and will be referred to as "raw". The final q^ 
spectrum is corrected for reconstruction effects by applying an unfolding algorithm to the measured 
(f' spectrum. The total B{B^ tt~£~^u) is given by the sum of the partial branching fractions 
AB{B^ — > 'K~ l'^u,q^). The q^ shape of the f^{q^) form factor is obtained from the normalized 
partial branching fractions I^B{q^)/B spectrum combined with two covariance matrices (one for the 
statistical errors and one for the systematic errors) which give the correlations among the values 
of /S.B{q^)/B measured in the different q^ bins. The measured /S.B[q^)/B spectrum is fitted to a 
model-dependent parametrization |S| of the f~^{q^) form factor. The model-independent /S.B{q^)/B 
spectrum and its correlation matrices are given explicitly to allow future studies with different 
f^i'f') parametrizations using the present data. The value of the CKM matrix element \Vub\ is 
then derived from the form factor calculations, combined with the measured lS.Bs. 

The main innovation of this analysis is the use of a loose neutrino reconstruction technique which 
yields a much higher signal reconstruction efficiency than in past measurements, while keeping the 
systematic errors at a relatively low level. The higher yield allows the utilization of a large number 
of (f' bins and the determination of the background composition using several independent fit 
parameters. In addition, we use g^-dependent cuts and estimate the f^{q^) shape systematic error. 

The data used in this analysis were collected with the BABAR detector at the PEP-II asymmetric 
e"'"e~ collider. The BABAR detector is described elsewhere jU]. The following samples are used: 206.4 
fb~^ integrated luminosity of data collected at the T(45') resonance, corresponding to 227.4 million 
BB decays; 27.0 fb~^ integrated luminosity of data collected approximately 40 MeV below the 

^Charge conjugation decays are implied throughout this paper. 
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T{4:S) resonance (denoted "off-resonance data" liereafter); standard BABAR Monte Carlo (MC) 
simulation using GEANT4 ^0] and EvtGen 1.64 million — > TT~i'^h' signal events using 
the FLATQ2 generator jl2j : 2.02 billion generic B^B^ and B+B' events, and 2.44 billion generic 
uu/dd/ss/cc and t^t~ "continuum" events. 

2 Analysis Method 

Values of \Vub\ have previously been extracted from B^ —>■ ir^t^v measurements by CLEO [TsllTl]. 
BABAR ^Sl El and BELLE . Our analysis is based on a novel technique denoted "loose neutrino 
reconstruction" . The main motivation for implementing this technique is to maximize the extracted 
B^ 7r~£^i' signal yields in order to measure the f~^{q^) shape as precisely as possible. 

Even though B mesons are always produced in pairs at the T(4S') resonance, a major feature 
of the neutrino reconstruction technique is that the decay of the non-signal B is not reconstructed. 
Instead, the signal B mesons are directly identified using the measured vr^ and £^ tracks together 
with the events' missing momentum as an approximation to the signal neutrino momentum |13[ I14 | 
[TB] . The neutrino four-momentum, Pmiss = {\Pmiss\-,'Pmiss)-, is inferred from the difference between 
the four- momentum of the colliding-beam particles, -Pbeams = (-Ebeams)Pbeams)i and the sum of the 
four-momenta of all charged and neutral particles detected in a single interaction, Ptot = {Etot-iPtot), 
such that Pmiss = Pbeams " Ptot- Compared with the tagged analyses described in Refs. ^| ITT] , 
the neutrino reconstruction approach yields a lower signal purity but a significant increase in 
the signal reconstruction efficiency. The new approach used in the loose neutrino reconstruction 
further increases this efficiency compared with the previous untagged analyses I14| II 5j by 
avoiding neutrino quality cuts (for example, a tight cut on the invariant missing mass to ensure the 
neutrino properties are well taken into account). Such cuts were required to allow the calculation 
of (f' = {Pi + Pjy)^. To obtain the values of q^, we use instead the neutrino-independent relation: 
(f' = {Pb — Ptt)"^ ■ Although this relation is strictly true and Lorentz invariant, it cannot be used 
directly because the value of Pb is not known. Only the value of Pt^, measured in the laboratory 
frame, and that of the T(4S') 4-momentum are known. Nevertheless, since the B momentum is 
small in the T{AS) frame, a common approximation is to boost the pion to the T(4S') frame and 
use the relation = {Pb — P-k)'^ in that frame, where the B meson is assumed to be at rest. 

However, a more accurate value of can be obtained in the so-called y-average frame 119j 
where the pseudo-particle Y has a 4-momentum defined by Py = {P-k + Pt)- The angle 9 by 
between the directions of the p*^ and py momenta^ in the T(4S') rest frame can be determined 
assuming energy-momentum conservation in a semileptonic B Yv decay. Its value is given 
by: cosOby = {'2,E*qEy — m\ ~ ^y) I ^\Pb\\Vy\)-' where mB,mY, E"^, Ey,Pb and py refer to the 
masses, energies and momenta of the B meson and the Y "particle". Thus, in the T(4S') frame, 
a cone is defined whose axis is given by the direction of the Y momentum with the half-angle 
subtended at the apex given by 6 by- The apex corresponds to the vertex formed by the Y and 
B momenta directions. The B momentum lies somewhere on the surface of the cone and thus 
its position is known only up to an azimuthal angle 4> defined with respect to the Y momentum. 
The value of q'^ in the Y-average frame is computed as follows: it first assumes that the B rest 
frame is located at an arbitrary angle (po, and the value of is calculated in that specific frame 
position. The values of , and q^ are then calculated with the B rest frame at cpi = (j)o + 90°, 
02 = 00 + 180° and 03 = 0o + 270°, respectively. The value of q^ in the Y-average frame is then 

®A11 variables denoted with an asterisk (e.g. p*) are given in the T(4S') rest frame; aU others are given in the 
laboratory frame. 
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defined as (f' = '^"^ ^'^ '^^ . Using more than four values of (pi does not significantly improve the 
resolution. 

The use of the y-average frame yields a resolution that is approximately 20% better than 
what is obtained in the usual T{4S) frame where the B meson is assumed to be at rest. We 
get an unbiased q'^ resolution of 0.52 GeV^/c^ when the selected pion candidate comes from a 

— > TT'l^v decay (Fig. which accounts for approximately 91% of our signal candidates after 
all the analysis selections. When a track from the non-signal B is wrongly selected as the signal 
pion, the resolution becomes very poor and biased. 
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Figure 1: resolution of B^ ir't'^v signal events obtained in the Y-average and T(45) frames 
after all analysis cuts and MC corrections. The very long tail arises when a track coming from the 
non-signal B is wrongly selected as the signal pion. The numbers of entries in the first and last 
bins correspond to the sum of all entries with Ag^ < -9.75 GeVVc"^ and Ag^ > 9.75 G^N'^/c^, 
respectively. 



We correct for our imperfect q^ resolution with a g^-unfolding algorithm. This algorithm was 
validated with statistically independent signal MC samples. After all selections, the total signal 
MC sample contains approximately 120000 events. Five thousand such signal events were used to 
produce the raw (f' and true q^ histograms. The remaining signal events were used to build the 
two unfolding matrices, using the simulated signal events reweighted [12] either to reproduce the 
f~^{q^) shape measured in Ref. or with the weights calculated in Ref. [7]. As illustrated in 
Fig. 12 the true and raw yield distributions differ considerably for various values of q^. However, 
the unfolded values of q^ match the true values within the statistical uncertainties of the unfolding 
procedure, independently of the signal generator used to compute the detector response matrix. 
This shows that the g^-unfolding procedure works as expected. 

To separate the B^ — > T^~i'^v signal from the backgrounds, we require two well reconstructed 
tracks that fulfill tight lepton and pion identification criteria. The electron (muon) tracks are 
required to have a momentum greater than 0.5 (1.0) GeV/c in the laboratory frame. We do 
not cut on the pion momentum because it is very strongly correlated with q^. The kinematic 
compatibility of the lepton and pion with a real B^ t:~1'^v decay is constrained by requiring 
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Figure 2: Validation of the g^-unfolding procedure. The true and raw yield distributions differ 
considerably for various values of q^. However, independently of the signal generator used to 
compute the detector response matrix, the unfolded values of match the true values within the 
statistical uncertainties of the unfolding procedure. 



that a geometrical vertex fit of the two tracks gives a probability greater than 0.01, and 
by requiring that — 1 < cos < 1. Note that cuts whose values depend on the measured value 
of (Fig. 3) give the best background rejection. Non-i?i? events are suppressed by several 
conditions: we require at least four charged tracks in each event; we require the ratio of the second 
to the zeroth Fox- Wolfram moments 21_ to be less than 0.5; we require the cosine of the angle 
between the y's thrust axis and the rest of the event's thrust axis, cos Othrust, to satisfy the relation^ 
cosethrust < 0.460+0.0576-^2 -0.00215-g^ (Fig. EJ; W6 rcQuirG the polar angle ctssocia-ted with Pmiss 
to satisfy the relation 2.7 rad > Omiss > (0.512 - 0.0162 • f + 0.000687 • g^) rad (Fig. EJ. Radiative 
Bhabha events are rejected using the criteria given in Ref. |22] and photon conversion events are 
vetoed. Finally, although the shapes of the q"^, AE and mss distributions in off-resonance data are 
very well reproduced by MC simulation in all lepton channels, there is an excess of nearly a factor 
of two in the yield values observed in data compared to the simulation in the electron/positron 
channels. We then require < 0.64 and > 0.35 for candidates in the electron and positron 
channels, respectively, where the z axis is given by the electron beam direction (Hj. This reduces the 
observed excess by removing additional radiative Bhabha events as well as "two-photon" processes 
which are not included in the simulated continuum. 

To reject background BB events, we require the Y candidates to have cos 6e < 0.85 and cos 6i > 
-0.938 + 0.0994 • q^ - 0.00384 • (Fig. EJ, where Oi is the hehcity angle of the W boson ^ 
reconstructed in the Y-average frame approximation. We reject J/V' — > IJ-'^ decays, which can 
often be mistaken for B^ — > 'k~ jl^v decays^, by removing candidates with 3.07 < my < 3.13 
GeV/c^. Of all the neutrino quality cuts utilized in Refs. j^lElEl) only the loose g^-dependent 

'^In the following relations, if is given in units of GeV^/c*. 

*This requirement is not necessary in the electron channel since the fake rate of charged pions by electrons is 
extremely low. 
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Figure 3: Cut functions used in the analysis, for the variables cos Othrust (a), cos 9i (b), Omiss (c) 
and M^i^3 (d). The sohd curves show the cut values as a function of q^. The histograms show the 
statistically optimal cut values obtained in 10 q'^ bins based on the minimization of the quantity 
y^{S + B)/ S in the AE-uies signal region, where 5* represents the simulated signal yield and B 
stands for the simulated background yield. 



criterion on the squared invariant mass of Pmiss is used: iVf^^^^ <Z (10.2 — 1.12 • q^ -t- 0.0625 * q^^ 
GeV^/c^ (Fig. El). We discriminate against the remaining backgrounds using the variables A.E = 

Pbeams 

s/2)/^/s and mss = \/('S/2 + PB ' Pbeams 
where ^/s is the total energy in the T{4S) center-of-mass frame. Only candidates with \AE\ < 
1.0 GeV and rriEs > 5.19 GeV/c^ are retained. When several candidates remain in an event after 
the above cuts, we select the candidate with cos 6*^ closest to zero and reject the others. This rejects 
30% of the combinatorial signal candidates while conserving 97% of the correct ones and reduces 
the sensitivity of our analysis to the simulation of the candidates' multiplicity. After all cuts, the 
total signal event reconstruction efficiency varies between 6.6% and 9.7%, depending on the bin, 
as shown in Fig. ^ 

To obtain the B^ tt^I'^u signal yield in each of the 12 reconstructed cp' bins, we perform a 2+1 
dimensional (AE'-ttt-eS) 'f') extended binned maximum likelihood fit based on a method developed 
by Barlow and Beeston [21]. The fitted data samples in each bin are divided into four categories: 
B^ — > signal and three backgrounds, b uiu, other BB, and continuum. The distinct 

structure of these four types of events in the 2-dimensional AE-rriES plane is illustrated in Fig. 
|SJ Since the correlation between AE and tties cannot be neglected and is difficult to parametrize, 
we use the AE-m-Es histograms obtained from the MC simulation as two-dimensional probability 
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Figure 4: Signal efficiency as a function of true q^. 



density functions (PDF). The simulated signal events are reweighted to reproduce the f~^{q^) 
shape measured in Ref. The (f' shape of the simulated non-BB continuum background is 

scaled to match the off-resonance data control sample containing both and /i^ events, while 
the scaling of the yields requires separate and /i^ samples. The fit of the MC PDFs to the 
experimental data gives the values of twenty parameters: twelve parameters for the twelve signal 
bins, three for the h uiv background, four for the other BB background, and one for the 
continuum background, as illustrated in Fig. El The number and type of fit parameters were chosen 
to provide a good balance between reliance on simulation predictions, complexity of the fit and 
total error size. The corresponding Ai? and mES fit projections in each bin for the experimental 
data are shown in Figs. [7| and |H1 We obtain a total signal yield of 5047 it 251 events, while for 
backgrounds the b — > uiv yield is 10015 it 548 events, the other BB yield is 32788 it 445 events, and 
the continuum yield is 9801 it 467 events. The fit has a value of 428/388 degrees of freedom. In 
the more restricted signal region (mss > 5.272 GeV/c^, jAE'l < 0.18 GeV), the total signal yield 
is 1340 lb 40 events and the total background yield 2527 it 55 events, for a signal/background ratio 
of 0.53 ±0.02. 

From the raw signal yields, the unfolded partial branching fractions AB{B^ — > vr^^^i^, g^) 
are calculated using the inverted detector response matrix given by the simulation and the signal 
efficiencies. The total branching fraction B{B^ — > 'K~t^v) is given by the sum of the partial 
/S.B{B^ — > 7r~£"'"z^, q^) thereby greatly reducing the sensitivity of the total branching fraction to the 
uncertainties of the f~^{q^) form factor values, which have a small but non-negligeable effect on the 
values of the efficiencies. 

To reduce the uncertainties in evaluating the f~^{q^) shape, instead of fitting the measured 
A;B(S° — >■ TT~t^i',q^) spectrum, we fit the normalized distribution, AB{q'^)/B, obtained by di- 
viding the measured AB spectrum by the measured value of the total branching fraction. With 
this approach, a number of correction factors cancel out, leading to a significant decrease in the 
systematic error. We fit the AB{q'^)/B spectrum using a PDF based on the f~^{q'^, a) parametriza- 
tion of Becirevic-Kaidalov jSj, which is proportional to the standard differential decay rate for a 
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Figure 5: Sum of AE-rriEs distributions over all 12 bins of for the four types of events used in 
the signal extraction fit after all cuts. 



semileptonic B decay to a pseudo-scalar meson (Xps): 

p. 2 ^ |p;M/+(g^«)P dTjB ^ Xps i^u,) 



where = W — ^ m^, and the f {q , a) function is: 

' " {l-qyml.yll-aqyml,) 

The value of /o cancels out in Eq. 1. Note that the data can also be used to extract the f~^{(f') 
shape parameter (s) using any theoretical parametrization, e.g. those of Refs. [SlUSlIini- The 
value minimized in the fit is defined in terms of the covariance matrix U to take into account the 
correlations between the measurements in the various q^ bins: 



Y: UB{qf)/B- Jj{q')dqA Ur/ (AB{q])/B - J J{q')dqA , (3) 
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Figure 6: binning used in the fit of the MC PDFs to the experimental data. 



where / 2 F{q^)dq^ denotes the integral of Eq. Qlover the range of the i*'^ bin and /„2 F{q^)dq^ = 
1. The central value of the parameter a, and its total error, are obtained using the total covariance 
matrix in Eq. 3. In the present case, in which the errors on AB{q'^)/B are more or less uniform 
across the q^ bins, using the statistical or the systematic covariance matrix in Eq. 3 yields the 
statistical or the systematic errors for a, respectively. Their quadratic sum is in fact consistent 
with the total error. The statistical covariance matrix is given directly by the fit to the signal. The 
systematic and total covariance matrices are obtained as described in the next section. 

3 Systematic Error Studies 

Numerous sources of systematic uncertainties have been considered. Their values are established by 
a procedure in which variables used in the analysis are varied within their allowed range, generally 
established in previous BABAR analyses. For the uncertainties due to the detector simulation, the 
variables are the tracking efficiency of all charged tracks (varied between it 0.7% and it 1.4%), the 
particle identification efficiencies of signal candidate tracks (varied between it 0.2% and it 2.2%), 
the calorimeter efficiency (used in the full-event reconstruction, and varied between it 0.7% and ± 
1.8% for photons, and up to it 25% for mesons) and the energy deposited in the calorimeter 
by mesons (varied up to it 15%). For the uncertainties due to the generator-level inputs to the 
simulation, the variables are the branching fractions of the background processes b uiv, b civ 
and D K^X as well as the branching fraction of the T(4S') B^B^ decay (all varied within 
their known errors 2^ except when the branching fractions have not been measured. In those cases, 
the branching fractions are varied by it 100% from their presumed central values). The B — > piv 
form factors are varied within bounds of it 10% at q^ = and it 16% at q^ax^ given by recent 
Light-Cone Sum Rules calculations ,27,, while the B D*iu form factors are varied within their 
measured uncertainties ^Hl; between ± 5.5% and ± 9.6%. To take into account an additional subtle 
effect on the uncertainty of the signal efficiency, the B^ — > 'K~t^v form factor shape parameter a 
is varied between its recently measured central value ^3] and that of the unquenched HPQCD 
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calculations a difference of 0.2. Finally, for the uncertainties due to the modelling of the 
continuum, there are variations in the continuum yields and in the g^, AE and rriEs shapes, as 
discussed in Section 2. 

The systematic errors are then given by the variation in the final values of the branching frac- 
tions when the data are re-analyzed with different values of the simulation parameters. For each 
source of uncertainty, we generate at least one hundred MC samples in which the simulation pa- 
rameters are varied according to a Gaussian standard deviation. This standard deviation is given 
by the range of the variations listed above. For each MC sample, the entire analysis is reproduced 
leading to new signal efficiencies, g^-unfolding matrices, AE-m^s PDFs and — > TT~i~^i' signal 
yields from a fit to the same data sample. The rms value of the resulting branching fraction dis- 
tribution is taken to be the value of the systematic error contributed by the source of uncertainty 
under study. The individual branching fractions are also used to generate two-dimensional AB{qf) 
versus AB{qj) distributions, for all (qf^qj) combinations. The linear correlation coefficient in each 
of these distributions is used to build the covariance matrix of the AB{q'^) measurements for each 
source of systematic error. The total systematic covariance matrix is then simply given by the sum 
of all the individual covariance matrices and is used to calculate the total systematic error on the 
branching fraction. The same procedure is repeated for the normalized branching fractions. The 
resulting total systematic covariance matrix yields in this case the total systematic error on the 
parameter a. All the statistical and systematic uncertainties are given in Tables A-1 and A-2 of 
Appendix A while the correlation matrices of the normalized branching fractions are presented in 
Tables A-3 and A-4. 



4 Results 

The values of the partial AB{B° tt-£+u, q^) and total B{B^ — > vr i^v) branching fractions 
are given in Table A-1, those of the normalized partial AB{q^)/B branching fractions are listed 
in Table A-2. In Table A-1, we also give the small uncertainties on the signal efficiency and q^- 
unfolding matrix due to the signal MC statistics. The total branching fraction error is due in 
large part to the photon and tracking efficiency systematic uncertainties. However, the use of the 
loose neutrino reconstruction did indeed reduce their impact |il5. . The systematic errors arising 
from the branching fractions and form factors of the backgrounds have been greatly reduced by 
the many-parameter fit to the background yields in the twelve bins of (f. As expected, the errors 
on AB{q'^)/B are mostly statistical. The value of the total B{B^ tt^I'^u) branching fraction 
obtained from the sum of the partial AB{B^ ■k~1~^v, q^) branching fractions is: 

B{B^ ^ TT-l+u) = (1.44 ± 0.08,tot ± Q.lOsyst) X 10-^ 

The normalized AB{q'^)/B distribution is displayed in Fig. [^together with the result of a 
f~^{q^) shape fit using the BK parametrization and theoretical predictions. We obtain a value of 
a = 0.53 lb 0.05 it 0.04. In Tabled we give the values and their associated probabilities for the 
four different calculations. These values were obtained by comparing, bin by bin, the data with 
the central values of the form-factor calculations, ignoring the theoretical errors. Our experimental 
data are clearly incompatible with the ISGW2 quark model. A more definitive choice among the 
remaining theoretical calculations must await a substantial increase in statistics. 

Monte Carlo studies have shown that there is no significant fit bias for the B{B^ 
but a small 3.8% bias in the f~^{<f') parameter a which has been incorporated in its systematic 
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Figure 9: Differential decay rate formula (Eq. ^ fitted to the normalized partial AB{q'^)/B spec- 
trum in 12 bins of q'^. The smaller error bars are statistical only while the larger error bars include 
statistical and systematic uncertainties. The BK parametrization (solid black curve) reproduces the 
data quite well (x^ = 8.8 for 11 degrees of freedom) with the parameter a = 0.53 ±0.05 ±0.04. The 
data are also compared to LCSR calculations W (dotted line), unquenched LQCD calculations [21 
(long dashed line), |4j (short dashed line) and the ISGW2 quark model (dash-dot line). 



error. Various cross-checks have also been performed. The results were obtained separately for the 
electron and muon decay channels, for the off-resonance data replacing the continuum PDF, for 
the different AE-niES and q'^ binnings, and for the variations of all the analysis cuts, one at a time. 
All the cross-check studies were found to be consistent with the final results. 

We extract \Vub\ from the partial branching fractions AB using \ Vub\ = \J ABj (r^A^), where 
= (1.536 ± 0.014) ps |21 is the lifetime and = Ar/|V^fep is the normalized partial decay 
rate predicted by various form factor calculations. We use the LCSR calculations for q^ < 16 
GeV^/c^ and the LQCD calculations for > iq GeV^c^- The results are shown in Table EJ The 
uncertainties of the form- factor normalization are taken from Refs. We obtain values 

of \ Vub\ ranging from 3.6 x 10"'^ to 4.1 x 10~^. For the most recently published unquenched LQCD 
calculation we obtain \Vub\ = (4.1 ±0.2stat ^ O-'^syst-oAFF) x 10~^. 
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Table 1: ^ values and associated probabilities for various QCD calculation predictions and for 
the ISGW2 model compared to our measured spectrum, for twelve degrees of freedom. 





stat error only 


stat+syst errors 


QCD calculation 




Proh{^) (%) 




Probix^) (%) 


ISGW2 


49.5 


< 0.01 


34.1 


0.07 


Ball-Zwicky jS] 


17.0 


14.9 


13.0 


37.2 


FNAL [11 


16.2 


18.2 


12.5 


41.0 


HPQCD ]3| 


14.2 


28.6 


10.2 


60.2 



Table 2: Values of \Vub\ derived from the form- factor calculations. The first two errors arise from 
the statistical and systematic uncertainties of the partial branching fractions. The third error comes 
from the uncertainties on due to the theoretical calculations. 





( GeV7c4) 




AB (10- 






AC 


{ps 






\Vub\ 


:io 






Ball-Zwicky Jj 


< 16 


1.07 




0.06 


± 


0.08 


5.44 


± 


1.43 


3.6 


± 


0.1 


± 


0.1 


+U.6 
-0.4 


HPQCD 0" 


> 16 


0.37 




0.04 


± 


0.03 


1.46 


± 


0.35 


4.1 


± 


0.2 




0.2 


+0.6 
-0.4 


FNAL m 


> 16 


0.37 




0.04 


± 


0.03 


1.83 


± 


0.50 


3.6 


± 


0.2 




0.2 


+0.6 
-0.4 


APE "F 


> 16 


0.37 




0.04 


± 


0.03 


1.80 


± 


0.86 


3.7 


± 


0.2 




0.2 


+1.4 
-0.7 



5 Summary 

The succesful development of the loose neutrino reconstruction technique shows that it is not always 
necessary to have very pure signal samples to control the systematic errors. The gain in statistical 
precision can overcome the negative features of large backgrounds. This technique could thus be 
used advantageously in future measurements, possibly those of other exclusive B — > X^ii^ decays. 

In the present analysis, we have obtained the total — > ir^i^v branching fraction from the 
values of the partial branching fractions measured in 12 bins of and the f~^{(f') shape parameter 
using the Becirevic-Kaidalov parametrization. We summarize these results in Table |31 together with 
the value of \Vub\ extracted from a recent calculation of the form factor jSj. 

Table 3: Summary of the main results. 



^(5" ^ TT-^+I/) = (1.44 ± Qmstat ± 0.10,^,4) X 10-4 

aBK = 0.53 ± 0.05,tet ± 0.04,j,,t 

iKfel = (4.1 ±0.2,f,f ±0.2,^,tl[i:^FF) X 10-3 



Our value for the B{B^ — > tt^I'^v) is the most precise measurement to date and, by itself, is 
of comparable precision to the current world average |28|: B{B^ —> Tr~£^i') = (1.34 it 0.08stat it 
O.OSsyst) X 10"^. The new value of the BK parameter is an improvement over our previous mea- 
surement a = 0.61 lb 0.09 15 (no systematic error quoted). 
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The errors in Table A-2 together with the correlation matrices of the statistical and systematic 
errors presented in Tables and lA-41 will allow the present data to be studied with different future 
f~^{q^) parametrizations. A simple calculation already shows that our data are incompatible 
with the predictions of the ISGW2 quark model. 
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Appendix A 

The values of the partial A.B{B^ t:~ i'^i>,q^) and total B{B^ -n^l^v) branching fractions 
are given in Table A-1, those of the normalized partial AB{q^)/B branching fractions are listed in 
Table A-2. All the statistical and systematic uncertainties as well as their associated correlation 
matrices are given in Tables A-1, A-2, A-3 and A-4. 
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Table A-1: Partial AB{B^ -K-l+v,q^) and total B{B^ TT'e+iy) (xlO'^) and their errors (xlO^) from all sources. 



iiitcr\als ( G(A^7'(:-^) 


0-2 


2-4 


4-G 


()-<S 


8-10 


10-12 


12-14 


14-lG 


1()-18 


18-20 


20-22 


22-26.4 


Total 


ry-<lG 


(r>iG 


fitted BF 


113.2 


151.1 


144.6 


131.8 


149.3 


115.7 


177.8 


87.7 


122.7 


63.3 


107.6 


78.7 


1443.6 


1071.2 


372.4 


fitted yield stat err 


22.9 


20.0 


19.9 


20.1 


22.8 


21.5 


24.3 


20.8 


23.4 


19.5 


20.0 


20.1 


83.4 


63.3 


44.1 


trie eS 


14.7 


1.5 


6.1 


3.3 


3.7 


3.4 


4.4 


4.0 


2.4 


4.6 


1.3 


4.4 


40.3 


39.7 


4.5 


7 cff 


15.2 


1.0 


5.3 


7.0 


3.3 


9.0 


7.3 


4.3 


3.2 


4.0 


2.8 


1.9 


56.9 


50.9 


6.7 




1.6 


1.0 


1.1 


1.4 


1.5 


0.8 


2.2 


1.6 


1.7 


1.5 


1.2 


1.8 


7.1 


5.2 


3.7 


Y PID & trk eff 


2.3 


2.9 


2.1 


2.8 


2.3 


1.7 


3.4 


1.5 


2.3 


0.9 


1.9 


1.8 


22.1 


15.8 


6.5 


continuum yield 


3.2 


0.6 


0.4 


0.1 


0.1 


0.2 


0.2 


0.6 


0.6 


1.3 


0.8 


2.1 


4.2 


2.2 


4.6 


continuum 


12.9 


2.3 


1.8 


1.2 


1.5 


1.3 


1.1 


2.0 


2.2 


3.6 


3.9 


8.4 


8.8 


8.8 


12.4 


continuum rriEs 


6.1 


0.6 


0.4 


0.1 


0.9 


0.6 


0.1 


0.7 


0.5 


1.1 


1.2 


2.0 


12.4 


7.7 


4.7 


continuum AE 


2.6 


2.2 


0.6 


1.4 


2.4 


0.5 


0.4 


1.0 


0.7 


1.4 


3.4 


3.3 


17.8 


9.2 


8.6 


K^BF 


8.1 


6.3 


8.2 


3.3 


3.9 


4.5 


3.4 


4.6 


5.5 


4.9 


4.3 


4.5 


51.6 


34.2 


17.5 


b BF 


3.0 


2.6 


1.5 


2.8 


6.2 


1.2 


4.3 


2.0 


1.7 


2.7 


1.5 


1.6 


17.0 


13.2 


5.0 


b ulu BF 


1.7 


1.8 


1.1 


1.0 


1.6 


1.6 


2.0 


1.5 


1.7 


2.9 


7.9 


7.5 


16.6 


9.6 


11.5 


T{AS) B"B" BF 


2.3 


3.3 


2.3 


2.0 


2.1 


1.9 


2.9 


1.1 


2.0 


0.8 


2.4 


1.4 


23.9 


17.7 


6.3 


B -> D*£iy FF 


1.7 


1.3 


0.2 


2.0 


4.0 


1.0 


2.6 


0.9 


1.0 


0.8 


0.8 


2.5 


12.5 


10.2 


2.7 


B p£u FF 


4.0 


1.2 


3.4 


1.7 


1.1 


1.7 


2.6 


3.9 


1.3 


1.9 


1.6 


3.5 


18.3 


14.4 


5.7 


^0 TT-e+v FF 


-1.2 


-0.0 


0.3 


0.1 


0.3 


0.0 


0.7 


-0.3 


-0.2 


-1.3 


1.6 


4.5 


4.7 


0.0 


4.7 


signal MC stat error 


1.8 


2.6 


2.4 


2.6 


2.4 


2.3 


2.6 


1.5 


1.7 


1.1 


1.3 


1.1 


5.5 


5.1 


2.2 


B counting 


1.2 


1.7 


1.6 


1.5 


1.6 


1.3 


2.0 


1.0 


1.4 


0.7 


1.2 


0.9 


15.9 


11.8 


4.1 


total syst error 


28.0 


9.8 


13.1 


10.6 


11.3 


11.7 


12.5 


9.7 


8.8 


10.3 


11.8 


15.6 


102.7 


83.0 


31.3 


total error 


36.2 


22.3 


23.8 


22.7 


25.5 


24.5 


27.3 


23.0 


25.1 


22.1 


23.3 


25.4 


132.3 


104.4 


54.1 



Table A-2: Normalized partial AB{q'^)/B (xlO^) and their errors (xlO^). 



intervals (GeVVc^) 


0-2 


2-4 


4-6 


6-8 


8-10 


10-12 


12-14 


1 4-16 


16-18 


18-20 


20-22 


22-26.4 


normalized partial BF 


78.4 


104.6 


100.2 


91.3 


103.4 


80.1 


123.2 


60.8 


85.0 


43.9 


74.5 


54.6 


fitted yield stat err 


1 Q 


io.y 


lo.o 




lo.o 


1 /I Q 


10. o 


1/1/1 


ID.Z 


lo.o 


lo.y 


1 /I n 

14. u 


sip'Tial MO fitat error 


1.3 


1.8 


1.6 


1.8 


1.7 


1.6 


1.8 


1.0 


1.2 


0.8 


0.9 


0.8 


track eff 


7.7 


2.3 


1.6 


0.8 


1.0 


0.6 


1.0 


1.2 


2.2 


2.1 


2.3 


4.3 


photon efF 


7.0 


3.4 


0.7 


2.1 


1.9 


3.0 


1.0 


0.8 


1.6 


1.2 


4.1 


1.8 


K°eSkE 


0.9 


0.9 


0.6 


0.9 


0.9 


0.4 


1.4 


1.1 


1.0 


1.0 


0.6 


1.2 


Y PID & trk eff 


2.3 


0.5 


0.2 


0.7 


0.3 


0.1 


0.7 


0.2 


0.4 


0.5 


0.2 


0.5 


continuum yield 


2.3 


0.3 


0.2 


0.2 


0.3 


0.1 


0.4 


0.3 


0.2 


0.8 


0.4 


1.4 


continuum 


8.9 


2.1 


1.6 


1.1 


0.6 


1.3 


0.7 


1.4 


1.1 


2.3 


2.3 


5.8 


continuum rriEs 


3.4 


0.4 


0.5 


0.7 


1.4 


0.2 


1.0 


0.2 


0.4 


0.4 


0.2 


0.9 


continuum AE 


0.9 


0.3 


0.8 


0.2 


0.5 


1.3 


1.8 


0.2 


0.7 


0.5 


1.5 


1.6 


B{D Kl) 


2.9 


3.3 


2.3 


1.9 


2.2 


2.6 


3.2 


2.5 


1.5 


2.3 


1.3 


1.3 


B{b civ) 


1.5 


2.8 


1.4 


1.3 


3.3 


1.0 


2.1 


1.3 


1.1 


1.6 


1.0 


0.8 


B{b-^ uiv) 


0.6 


1.3 


0.6 


0.9 


1.2 


0.7 


1.7 


0.6 


1.6 


2.0 


4.8 


4.9 


B{Y{'iS) B°BO) 


0.3 


0.6 


0.2 


0.2 


0.4 


0.1 


0.3 


0.3 


0.1 


0.6 


0.4 


0.2 


B D*iv FF 


0.6 


1.7 


0.8 


0.7 


2.0 


0.5 


0.8 


0.6 


1.3 


0.7 


0.3 


1.3 


B plv FF 


1.8 


0.6 


1.4 


0.6 


0.7 


0.3 


3.2 


2.0 


0.7 


1.4 


1.5 


1.8 


B^ TT-e+v FF 


1.1 


0.3 


0.1 


0.2 


0.1 


0.2 


0.1 


0.4 


0.4 


1.0 


0.9 


2.9 


total syst error 


15.1 


7.1 


4.5 


4.3 


5.7 


4.9 


6.4 


4.4 


4.5 


5.5 


7.8 


10.1 


total error 


21.9 


15.6 


14.5 


14.6 


16.8 


15.7 


18.0 


15.1 


16.9 


14.6 


15.9 


17.2 



Table A-3: Correlation matrix of the normalized partial /\B{(f')/B statistical errors. 



2-4. 1 

q intervals 

(^vjev /( ) 


0-2 


2-4 


4-6 


6-8 


8-10 


10-12 


12-14 


14-16 


16-18 


18-20 


20-22 


22-26.4 


0-2 


1.00 


-0.20 


0.12 


-0.00 


-0.01 


0.04 


0.04 


-0.01 


-0.00 


0.01 


0.01 


0.01 


2-4 


-0.20 


1.00 


-0.32 


0.14 


0.03 


0.01 


0.02 


-0.01 


-0.00 


-0.00 


0.00 


0.00 


4-6 


0.12 


-0.32 


1.00 


-0.31 


0.20 


0.05 


0.13 


-0.02 


-0.00 


-0.00 


0.00 


0.00 


6-8 


-0.00 


0.14 


-0.31 


1.00 


-0.22 


0.14 


0.08 


-0.01 


-0.00 


-0.00 


-0.00 


-0.00 


8-10 


-0.01 


0.03 


0.20 


-0.22 


1.00 


-0.23 


0.19 


-0.03 


0.00 


-0.00 


-0.00 


-0.01 


10-12 


0.04 


0.01 


0.05 


0.14 


-0.23 


1.00 


-0.02 


0.02 


-0.00 


0.00 


-0.00 


0.00 


12-14 


0.04 


0.02 


0.13 


0.08 


0.19 


-0.02 


1.00 


-0.24 


0.01 


-0.02 


-0.00 


-0.00 


14-16 


-0.01 


-0.01 


-0.02 


-0.01 


-0.03 


0.02 


-0.24 


1.00 


0.00 


0.11 


-0.03 


-0.00 


16-18 


-0.00 


-0.00 


-0.00 


-0.00 


0.00 


-0.00 


0.01 


0.00 


1.00 


0.01 


-0.04 


-0.01 


18-20 


0.01 


-0.00 


-0.00 


-0.00 


-0.00 


0.00 


-0.02 


0.11 


0.01 


1.00 


-0.18 


-0.11 


20-22 


0.01 


0.00 


0.00 


-0.00 


-0.00 


-0.00 


-0.00 


-0.03 


-0.04 


-0.18 


1.00 


-0.01 


22-26.4 


0.01 


0.00 


0.00 


-0.00 


-0.01 


0.00 


-0.00 


-0.00 


-0.01 


-0.11 


-0.01 


1.00 



Table A-4: Correlation matrix of the normalized partial AB{q^)/B systematic errors. 



intervals 
(GeV7c4) 


0-2 


2-4 


4-6 


6-8 


8-10 


10-12 


12-14 


14-16 


16-18 


18-20 


20-22 


22-26.4 


0-2 


1.00 


-0.58 


0.04 


-0.20 


-0.27 


0.19 


-0.21 


0.12 


-0.48 


0.15 


-0.40 


-0.48 


2-4 


-0.58 


1.00 


0.09 


-0.16 


-0.01 


-0.35 


0.12 


-0.16 


0.38 


-0.39 


0.34 


0.16 


4-6 


0.04 


0.09 


1.00 


-0.02 


-0.26 


0.16 


-0.25 


0.12 


0.06 


-0.05 


-0.28 


-0.15 


6-8 


-0.20 


-0.16 


-0.02 


1.00 


0.28 


0.46 


0.33 


-0.05 


-0.03 


-0.20 


-0.37 


-0.11 


8-10 


-0.27 


-0.01 


-0.26 


0.28 


1.00 


-0.38 


0.52 


-0.35 


0.01 


-0.33 


0.15 


-0.06 


10-12 


0.19 


-0.35 


0.16 


0.46 


-0.38 


1.00 


0.03 


0.33 


-0.22 


0.14 


-0.53 


-0.28 


12-14 


-0.21 


0.12 


-0.25 


0.33 


0.52 


0.03 


1.00 


-0.53 


0.12 


-0.34 


0.03 


-0.35 


14-16 


0.12 


-0.16 


0.12 


-0.05 


-0.35 


0.33 


-0.53 


1.00 


-0.33 


0.36 


-0.17 


-0.10 


16-18 


-0.48 


0.38 


0.06 


-0.03 


0.01 


-0.22 


0.12 


-0.33 


1.00 


0.07 


0.12 


0.04 


18-20 


0.15 


-0.39 


-0.05 


-0.20 


-0.33 


0.14 


-0.34 


0.36 


0.07 


1.00 


-0.29 


-0.04 


20-22 


-0.40 


0.34 


-0.28 


-0.37 


0.15 


-0.53 


0.03 


-0.17 


0.12 


-0.29 


1.00 


0.17 


22-26.4 


-0.48 


0.16 


-0.15 


-0.11 


-0.06 


-0.28 


-0.35 


-0.10 


0.04 


-0.04 


0.17 


1.00 
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